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ABSTRACT 

PSR B0628-28 is a radio pulsar which was first detected in the X-ray band by 
ROSAT and then later observed with Chandra and XMM-Newton. The Chan- 
dra observation yielded an X-ray luminosity two orders of magnitude higher than 
what is expected for spin-powered pulsars, also there were no pulsations detected. 
The XMM- Newton observation, however, reveals pulsations at the expected radio 
period, P = 1.244 s. The simultaneously analyzed spectra also give a luminos- 
ity (in cgs) logL x = 30.34, which is ~ 350 times greater than what would be 
expected from the correlation between L x -E. 

Subject headings: pulsars:individual (PSR B0628-28) - - stars:neutron - - X- 
rays: stars 

1. Introduction 

PSR B0628-28 is an old pulsar that has a characteristic age (r = P/2P) of 2.8 Myr. 
Since older pulsars have radiated away their initial heat content and have relatively low 
rotational energy loss (E = Iflfl) , one would expect to detect X-ray radiation coming either 
from the reheating of the surface (e.g. see Tsuruta 1998; Shibazaki & Lamb 1989) or by 
heated polar caps (Halpern & Ruderman 1993). In young pulsars, however, thermal radiation 
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from the neutron star surface is dominated by the nonthermal component from the neutron 
star magnetosphere, whose spectrum can be described primarily by a power-law model. We, 
also, know from X-ray observations that the pulsar's nonthermal luminosity (Lx) shows a 
correlation with the spin-down power E. For instance, Possenti et al. (2002) found a best fit 
based on 37 pulsars. Although the variance is large, pulsars follow a general trend that can be 
formulated as logL^f 10 = 1. 34 log E — 15.34. This set included only three old (r > 10 6 years) 
pulsars. The reason for having only three is two fold. The first one is that older pulsars 
are less active and relatively cooler, hence most of them were not detected until Chandra 
and XMM- Newton. The detected ones have low countrates, consequently statistically poor 
spectra. This makes it challenging to distinguish between spectral models. The other reason 
is that it is expected that surface emission from old neutron stars will be dominant, resulting 
in spectra to be more like a blackbody. Hence, the X-ray radiation from these pulsars are 
not expected to be characterized by a power-law model. However, recent observation of 
PSR B0943+10 (Zhang et al. 2005), PSR B0823+26, PSR B0950+08, PSR J2043+2740 with 
XMM-iVewton (Becker et al. 2004) and of PSR B2224+65 with Chandra (Zavlin & Pavlov 
2004) has unfolded a new perspective in understanding old neutron stars. Not only many 
previously undetected sources were detected but also it has been possible to distinguish 
between different spectral models for some of these old pulsars. These old pulsars have 
spectra harder than expected and is best described by a power-law and/or a polar cap 
(black-body) model, as opposed to a black-body model representing the thermal emission 
from the whole neutron star surface. The integrated model flux of these sources convert to 
unusually high X-ray efficiencies (L x /E). 

PSR B0628-28 is the longest period radio pulsar detected in X-rays. However, physical 
parameters inferred from radio observations are not in any means extreme. The Chandra 
data yielded a luminosity two orders of magnitude greater than what is expected for a spin- 
powered pulsar (Ogelman & Tepedelenlioglu 2004). The varying efficiency of converting 
spin-power in to X-ray luminosity in pulsars can be explained by geometrical effects. But, 
non of these can amount to such a great excess flux in the X-ray band. In this paper we 
try to address the reasons for this extreme luminosity. We describe the XMM- Newton and 
Chandra data in § 2 and present the result of our timing and spectral analysis in § 3. We 
discuss the implications of the analysis results in § 4. 

2. Chandra and XMM-IVewton observations of PSR B0628-28 

2.1. Chandra 



PSR B0628-28 was observed twice, first on 2001 November 04 and again on 2002 March 
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25 for 2000 s and 17000 s, respectively. For both observations photons were collected using the 
Advanced CCD Imaging Spectrometer (ACIS). Data were collected in the nominal timing 
mode, with 1.141 s exposures between CCD readouts. We reprocessed the Level 1 event 
data to correct the detrimental effects of charge transfer efficiency. The imaging, timing 
and spectral analysis presented here were done only on the 17 ks observation, the 2.0 ks 
observation was disregarded due to an interruption by a large solar storm. The background 
countrate during this solar storm increased by a factor of ~20. The net source countrate 
during the first observation period was 0.011 ±0.001 counts per second (cps) where the error 
is in the 68% confidence range. Since by taking into account the 2.0 ks observation we gain 
only ~27 counts, which is not enough to improve our statistics significantly, we preferred to 
disregard these counts, which potentially can be misleading. 

The measured PSF of PSR B0628-28 is consistent with the ACIS point-source response, 
hence the ACIS image reveals a point-like X-ray source at the pulsar position. The Chandra 
position of PSR B0628-28 is a = 06 h 30 m 49'.'43, 5 = -28°34'43'.'60 (J2000.0), which con- 
sidering 0'.'5 rms error and the ~0'.'6 absolute astrometric accuracy of Chandra, is in good 
agreement with the radio position a=06 h 30 m 49"53, 5 = -28°34'43'. / 60 (J2000.0), which was 
taken from the ATNF pulsar catalogue 1 . 



2.2. XMM.-Newton 

PSR B0628-28 was observed with XMM-Newton on 2004 February 28 for a total ontime 
of 48 ks. MOS1/2 were both operated in imaging (PrimeFullWindow) mode and the medium 
filter was used. During the EPIC-pn exposure the thin filter was used and the detector was 
operated in imaging (PrimeLargeWindow) mode for 47 ks. The temporal resolution achieved 
with this choice of science modes were 2.6 s and 43 ms for MOS1/2 and pn, respectively. 

The background of EPIC camera is known to be effected by soft proton flares. In order 
to screen for times of high background we inspected the light curves of MOS1/2 and pn data 
separately at energies above 10 keV. Both sets of data were effected by high background. The 
light curves were formed with 100 s bins and after inspection we rejected bins with countrate 
greater than 0.4 cps. The removal of high background time intervals from the data leaves 
us with effective exposure times of 42.5 ks and 33.3 ks for MOS1/2 and pn, respectively. 
The pulsar is clearly detected in the EPIC image at a = 06 h 30 m 49'.'48, 5 = -28°34'43"10 
(J2000.0) which differs from the radio position by only 0'.'8 well within the 2"-3" uncertainty of 
the EPIC absolute astrometry. The shape of the radial profile of the source is also consistent 
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with that expected for a point-like source. 

3. Results 

3.1. Timing 

For searching pulsations from PSR B0628-28 Chandra ACIS and XMM-Newton EPIC- 
MOS data were not suitable due to their limited temporal resolution. The sampling frequency 
in both cases set by the detector read out rate give a Nyquist frequency greater than the 
pulsar frequency (see Table 1). Thus we used EPIC-pn data that has 43 ms timing resolution. 

We extracted source plus background photons from a 30" radius circle centered at the 
pulsar position, which encircles about 85% of all detected source counts. The extraction 
region contained 1047 counts of which 16% is background. The photon arrival times were 
solar system barycenter corrected. The pulsar's spin parameters (Table 1) are well known 
from radio observations and can be extrapolated to the mean epoch of the XMM-Newton 
observation: MJD=53063.339. Around this predicted pulsar frequency, we then generated a 
periodogram using the Zf-statistic. The X-ray periodogram is shown in Figure 1. We found 
a peak at / = 0.80358444 ± 0.00000112 Hz which is consistent with the extrapolated pulsar 
frequency / = 0.80358551 Hz. The Zf for this peak is 40.4, which has a probability of chance 
occurrence of 1.69 xl0~ 9 . The pulse profile of PSR B0628-28 over the whole energy band 
(bottom panel in Figure 2) is broad and single-peaked with a pulse fraction of f p = 35 ±12%. 
Here we defined the pulse fraction as (Cmax _ C m i n )/(Cmax + ^min)' where Cmax and 
C m j n are the maximum and minimum counts per bin, respectively. 

We also looked at the pulse profile of the pulsar at different energy bands (upper panels 
of Figure 2). Most of the counts are in the soft band, 67% and 87% of all counts are in 
the 0.2-1 and 0.2-2 keV energy bands, respectively. Pulsed fraction in each selected energy 
band seems to be consistent with each other within the la errors associated. The pulse 
shape does not seem to be energy dependent which would suggest that all pulsed X-rays are 
coming from the same region. 

3.2. Spectral 

The pulsar's energy spectrum was extracted from the MOS1/2 data by selecting all 
events detected in a circle of radius 50" centered on the pulsar position. This region includes 
90% of all event from the pulsar. Due to a source located close to the pulsar we extracted the 
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background spectrum from a nearby circular region with radius 87". For the EPIC-pn data 
we extracted the spectrum from a circle centered on the pulsar with radius 30" (includes 
80% of source counts). PSR B0628-28 was being located very close to the chip boundary 
precluded the extraction of the background spectrum from an annular region around the 
pulsar. Hence, we used an off-source circular region with radius 66". 

To extract the spectrum from the Chandra data we used a circular region centered on 
the pulsar position with 2" radius. This region contains 95% of all source counts. The 
background spectrum was extracted from an annulus of radii 3"< r <50". 

In total, the extracted spectra include 780 EPIC-pn source counts and 754 EPIC- 
MOS1/2 source counts. Both spectra were binned so that each bin contained minimum 
25 counts per bin. Chandra data had a total of 184 source counts. The extracted pho- 
tons were binned and regrouped such that each fitted spectral bin contained minimum of 20 
counts. All three extracted spectra were then simultaneously fitted with model spectra. 

Among the single-component spectral models, an absorbed power-law model gave the 
statistically best representation (x 2 = 56.2 for 62 degrees of freedom (dof)) of the observed 
spectrum. A single black-body (x 2 = 93 for 62 dof) did not give a statistically acceptable 
fit. This fit when absorbing column is left to vary also yields a very low iVjj. The best fit 
power-law spectrum and residuals are shown in Figure 3. 

The power-law model yields a column density of iVjj = 1.381q 23 x 10 21 cm -2 , a photon 
index T = 3.20±g;||, and a normalization of 1.73l£^ x 10~ 5 photons cm" 2 s" 1 keV" 1 at E = 1 
keV. The errors are the upper and lower bounds of the la confidence range. The normal- 
ization converts to an unabsorbed energy flux of fx~ 10 — 8-6lla33 x 10~ 15 erg cm' 2 s _1 in 
the 2-10 keV band. Given the distance of d — 1.45 kpc this yields an X-ray luminosity of 
L|" 10 = 2.171° x 10 30 erg s 1 . This luminosity implies a rotational energy to X-ray energy 
conversion factor L x / E = 0.015 within 2-10 keV band. 

It is natural to assume that, in addition to the magnetospheric emission, thermal emis- 
sion from polar caps or from the surface due to reheating of the superfluid interior contributes 
to the observed X-ray flux. In order to explore this possibility and how it represents the 
data we used a two component model, thermal and magnetospheric, to fit the spectra. As 
a first approach we let every parameter vary. The fit gives a blackbody temperature of 
T = 3.28lgg2 x 10 6 K and effective radius R = 591^ m. The power-law photon index does 
not change significantly and has a value of T = 2.98lo'g5. The bolometric luminosity of the 
thermal component is L-u j = 2.87 x 10 30 erg s _1 whereas the total luminosity in the 2-10 
keV is Llr w = 1.67±g;§2 x 10 30 erg s" 1 . The nonthermal X-ray luminosity in the 2-10 keV 
band converts to an X-ray efficiency of L x /E = 0.01. The hydrogen column density is not 



- 6- 



well bound N E = 0.62+°;^ x 10 cm but still is consistent with the estimate obtained 
from the single power-law fit. The quality of the composite model (x 2 = 54.5 for 62 dof) 
is slightly better than that obtained from single power-law model. However, both fits are 
statistically acceptable and addition of a new model will always tend to decrease the x 2 ■ 

4. Discussion 

With the observations of new sensitive X-ray telescopes like XMM-Newton and Chan- 
dra we are uncovering a class of new sources; X-ray luminous old-radio pulsars. Up to 
now there are only seven such sources that have been detected. Namely, these are PSRs 
B2224+65, J2043+2740, B1929+10, B0823+26, B0950+08, B0943+10 and B0628-28 (see 
Zavlin & Pavlov 2004; Becker et al. 2004; Zhang et al. 2005; Ogelman & Tepedelenlioglu 
2004). Of these six sources only two of them (PSRs B0950+08, B0628-28) have high quality 
spectra so that one can distinguish between spectral models (e.g. thermal vs. nonthermal). 
Becker et al. (2004) suggested that the three pulsars observed with XMM-Newton (PSRs 
B0950+08, 0823+26, J2043+2740) all have spectra described better by a single power-law 
model, indicating nonthermal emission is dominant. For these pulsars, however, there are al- 
ternate representations that do give statistically and physically acceptable results (Zavlin & 
Pavlov 2004). For example, Zavlin & Pavlov (2004) argues that the soft part of the spectrum 
of PSR B0950+08 can be interpreted as radiation from heated polar caps on the neutron 
star surface covered with a hydrogen atmosphere. For PSR B2224+65 Zavlin & Pavlov 
(2004) arrives at the conclusion, from the analysis of the Chandra ACIS observation, that 
the pulsar definitely shows non-thermal emission. All these pulsars have X-ray luminosities 
greater than what is predicted from the correlation between TJ^^-E (Possenti et al. 2002). 
For example, in the case of PSR B0628-28 and PSR B2224+65 the luminosities inferred from 
spectral fits are greater then the so-called "maximum efficiency line" derived by Possenti et 
al. (2002) such that all pulsar lie below this line in the logL x -log-E plane. PSR B0628-28 
has been known to be an exceptional emitter and with the detection of pulsations from this 
source at the radio frequency has left no doubt that it is the X-ray counterpart of the pulsar. 
These luminous old pulsars (in particular PSR B0628-28) seem not to follow the trend their 
possible progenitors do. Which suggests that pulsars become more X-ray efficient as they 
grow older, given that r is inversely proportional to E. From multi-wavelength observations 
Zharikov et al. (2004) have arrived at the same conclusion. However, Zavlin & Pavlov (2004) 
suggest that the most plausible reason for the observed high nonthermal X-ray efficiencies 
associated with the old pulsars is the geometrical effects. We discuss these below. 

As mentioned in § 1, Possenti et al. (2002) used data from 39 X-ray emitting pulsars 
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to find a best fit that describes the correlation between logL x -log-E. These 39 sources were 
divided in to subcategories as; millisecond pulsars (P <10 ms), Crab-like (r ~ 10 4 years), 
Vela-like (r ~ 10 4 -10 5 years), Geminga-like (r > 10 5 years, where substantial amount of 
the X-ray flux comes from the internal cooling) and finally Old-pulsars (r > 10 6 years). In 
Figure 4 we show a similar plot. In this plot we only include Crab, Vela and Geminga- 
like pulsars. We did not include millisecond pulsars because they are re-cycled and do not 
represent the naturally aging pulsars with the same E. Also when converting countrates 
in to luminosities, Possenti et al. (2002) adopted a power-law spectrum with T = 2 for all 
millisecond pulsars. Hence, the luminosities do not come from spectral analysis and rather 
from assumptions based on other pulsars. Also we excluded the three Old-pulsars because 
their fluxes were obtained by scaling their ROSAT countrates to that of PSR B1929+10. 

This newly formed set of pulsars represent the evolution of pulsars on the log L x -log E 
plane. Where older pulsars are on the bottom left (low E) and young Crab-like pulsars are on 
the top right (high E). Using only these three subclasses (26 sources) we performed a linear 
fit of the form log/oc = a\ogE + b (see Figure 4). Due to exclusion of the stated pulsars 
our fit yields a line with a steeper slope (a ~ 1.5). We also identified the line of maximum 
efficiency as the line for which every pulsar lies underneath. We then overlayed only two of 
the seven old pulsars, with well known spectra. We should note that, however, with their 
luminosities derived from the power-law fits, PSR B0943+10 (Zhang et al. 2005) and PSR 
B2224+65 (Zavlin & Pavlov 2004), have efficiencies exceeding the maximum efficiency 

From Figure 4 it is apparent that old pulsars have very high nonthermal X-ray efficien- 
cies. This is contrary to what Possenti et al. (2002) has found and to what has been found 
when millisecond pulsars were excluded. The observed luminosities for younger pulsars also 
show large deviations, from this dependence. But in their case the deviation is symmetric 
around the trend line. The scatter that young pulsars exhibit could be due to uncertainties 
in pulsar distances, spread in the orientations of magnetic and rotational axes versus the line 
of sight. For example, seeing a certain fraction of the beam would result in lower inferred 
X-ray efficiency and vice versa. 

Although there should be a thermal contribution to the overall luminosity, this effect 
will not change the nonthermal luminosity significantly. For instance, when we subtract the 
thermal luminosity, obtained for PSR B0628-28 (see § 3.2), from the total, the nonthermal 
luminosity only changes by a factor of 0.5. The resulting luminosity is still too big and gives 
a nonthermal X-ray efficiency 2-3 orders of magnitudes greater than the maximum efficiency. 

The quality of the observed spectrum for these old pulsars are in general not high enough 
to distinguish between a thermal and a non-thermal model, or a combination of both. The 
current set of old pulsars all either show a non-thermal emission or the statistics is not high 
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enough to rule it out. For example, Becker et al. (2004) suggests that the X-ray emission from 
old, nonrecycled rotation powered pulsars is dominated by non-thermal radiation. When the 
data is fit by a power-law model, which in all cases is statistically acceptable, the obtained 
power-law indices are on average larger than the younger rotation powered pulsars. The 
mean of the photon indices of the 26 pulsars (Figure 4) is ~1.9 as opposed to ~2.5 for the 
old pulsars. The power-law fits suggests that the spectra of these old pulsars being steeper 
should result in lower luminosities in the 2-10 keV band. However, we observe an inverse 
effect where not only the spectra are steeper but also they have higher luminosities. Further 
X-ray observations of old radio pulsars should help us understand these puzzling features. 
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Table 1. Properties of PSR B0628-28 derived from radio observations 



Parameter Value 



Frequency (Hz) 0.80358811986 

Frequency derivative (Hz s _1 ) -4.59962xl0 _1 

Epoch (MJD) 46603.0 

Spin-down age (10 6 yr) 2.77 

Spin-down energy (10 32 ergs s -1 )... 1.5 

Inferred Magnetic Field (10 12 G)... 3.0 

Dispersion Measure (pc cm~ 3 ) 34.5 

Distance (kpc) 1.45 



Note. — Distance is inferred from dispersion-measure 
(Cordes & Lazio 2002). 



-10- 




Frequency (Hz) 



Fig. 1. — The power spectrum of the timing data around the expected frequency. The inset is the 
expanded view ±500 fiRz around the peak, which is located at / = 0.80358551 Hz. 
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Phase 



Fig. 2. — X-ray light curves of PSR B0628-28 extracted from the EPIC-pn data in four energy 
bands, with the values of the intrinsic pulsed fraction f p and its la errors. Two phase cycles are 
shown for clarity. For the expression used for the pulse fraction see text. 
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Energy (keV) 

Fig. 3. — Energy spectrum of PSR B0628-28 as observed with pn (upper spectrum), M0S1/2 
(lower spectrum) and ACIS-S (middle spectrum). The lines are the best fit power-law model. The 
lower panel is fit residuals. 
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Fig. 4. — The X-ray luminosity in the band 2-10 keV versus the spin-down power for the 26 sources 
together with two old pulsars. The value for the X-ray luminosity of PSR B0950+08 is taken from 
Becker et al. (2004) and the rest are taken from Possenti et al. (2002). The solid line is the best fit 
to the set of 26 pulsars. Other two lines are labeled in the legend. 



